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Abstract 

Background: Ocular melanosis of Cairn terrier dogs is an inherited defect characterized by progressive 
pigmentation of both eyes which can result in glaucoma and blindness. Pedigree analysis suggests the trait has an 
autosomal dominant mode of inheritance. We selected 1 1 potential candidate genes and used an exclusion 
analysis approach to investigate the likelihood that one of the candidate gene loci contained the Cairn 
terrier-ocular melanosis locus. 

Results: Two polymorphic loci were identified within or close to each candidate gene. Genotyping of at least 10 
ocular melanosis Cairn terriers for each marker showed that there was no single shared allele for either of the 
two polymorphic markers identified in ASIP, COMT, GPNMB, GSK3B, LYST, MC1R, MITF, SILV, TYR, TYRP1, and TYRP2. 
This is strong evidence to exclude each locus as the site of the ocular melanosis mutation (probability of a false 
exclusion calculated for each gene ranged from 1.59 x 10" 4 to 1 x 10" 9 ). 

Conclusions: None of the 1 1 potential candidate genes selected are likely to be the gene locus for ocular 
melanosis in Cairn terriers. 
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Background 

Ocular melanosis (OM) is a hereditary condition in the 
Cairn terrier breed of dog that is characterized by prolif- 
eration and migration of pigmented cells starting in the 
anterior uvea of both eyes. The pigmented cells result in 
a diffuse expansion of the root of the iris and are also 
shed into the aqueous humor. They accumulate in the 
anterior sclera which develops pigment plaques that pro- 
gressively increase in size, possibly reaching that site via 
the conventional aqueous drainage pathways [1,2]. Even- 
tually the aqueous drainage pathways become blocked 
with pigment results in secondary glaucoma and vision 
loss. The pigmented cells appear to be primarily melano- 
cytes although some admixture with melanophages 
(macrophages that have ingested melanosomes) does 
occur [3]. The pigmented cells also reach the posterior 
segment of the eye and over time can be seen to slowly 
migrate into and obscure the tapetum [2]. The posterior 
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sclera also becomes pigmented and pigment can be 
found in the meninges around the optic nerve. Trans- 
formation to develop solid melanocytic tumors occurs in 
some cases and was reported in 3 dogs in a series of 114 
ocular melanosis affected Cairn terriers [3]. 

Cairn terrier OM shares some similarities with human 
pigment dispersion syndrome (PDS). PDS patients accu- 
mulate pigment within the trabecular meshwork, and a 
subset (approximately 35% by 35 years following PDS 
diagnosis) develop pigmentary glaucoma (PG) [4], This 
condition accounts for 1-1.5% of human glaucoma cases 
seen in the Western world [5]. Although the PDS causa- 
tive mutation/s has yet to be identified, the phenotype 
has been linked to the telomere region (7q35-q36) of 
human chromosome 7 [6]. 

Rodent models with pigment changes similar to human 
PDS and PG, and pseudoexfoliation have been described. 
In the DBA/2 J mouse, mutations in GPNMB and TYRP1 
result in iris pigment dispersion and iris stromal atrophy 
respectively and mice homozygous for recessive alleles 
at both loci have earlier onset and more severe disease 
[7,8]. Clinically this mouse model closely resembles 
Cairn terrier OM. More recently a mutation in LYST 
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has been shown to cause iris transillumination defects 
in mice [9]. 

OM in Cairn terriers represents a potentially im- 
portant large animal model for investigation of human 
PDS and PG. Pedigree analysis indicates that Cairn 
terrier OM is likely to be inherited as an autosomal 
dominant trait [2]. We tested the hypothesis that one 
of 11 candidate genes has suffered a mutation that 
causes OM phenotype by examination of allele sharing 
of genetic markers close to or within each gene. Using 
an association-based form of exclusion analysis (as op- 
posed to pedigree-based exclusion analysis), we pro- 
vide evidence that strongly supports the exclusion of 
these 11 candidate genes from association with OM in 
the Cairn terrier. 



Results 

Candidate genes 

Eleven potential candidate genes for ocular melanosis 
were selected: agouti signaling protein (ASIP), catechol-o- 
methyltransferase (COMT), glycoprotein NMB (GPNMB), 
glycogen synthase kinase 3-beta (GSK3B), lysosomal traf- 
ficking regulator (LYST), melanocortin 1 receptor (MCIR), 
microphthalmia transcription factor (MITF), silver (SILV), 
tyrosinase (TYR), tyrosinase related protein 1 (TYRP1), and 
Tyrosinase related protein 2 (TYRP2). For each gene two 
polymorphic markers were identified (11 single nucleotide 
polymorphisms (SNPs), 9 microsatellites (MS), 2 insertion/ 
deletions (in/dels)) that could be used to test for associ- 
ation of the candidate locus with Cairn terrier OM. Details 
of the candidate genes and markers are shown in Table 1. 



Table 1 Candidate genes for Cairn terrier ocular melanosis with markers, chromosomal positions, and PCR primers 



Gene, Protein, Position 3 


Marker 13 [Restriction enzyme] 


Marker Position 0 


Distance (kb) d 


P values 8 


LYST Lysosomal Trafficking Regulator 


Microsat 1 (TTTTC) 1 80TTTC) 1 1 


1 § 


345 


1.00 x 10" 5 


chr4:7,1 28,220-7,294,031 


Microsat 2 (GAAA)5 (G AAA) 1 3 


2 § 


368 




MCIR Melanocortin 1 receptor 


SNP 1 (BICF2P987741) 


2 A 


560 


1.34X 10~ 5 


chr5:66,692,398-66,693,344 


SNP 2 (BICF2S23213233) [Ssil] 


1 A 


24 




SILV Silver 


Microsat 1 (GAAA)7(GAAA)1 7 


2 § 


60 


1.00 x 10" 5 


chrl 0:3,273,996-3,279,352 


Microsat 2 (mC)15(TTCC)13 


1 § 


226 




TYRP1 Tyrosinase related protein 1 


SNP 1 (BICF2S23051528) [BC1I] 


1 A 


224 


2.15 x 10" 5 


chrl 1:36,344,71 2-36,361,793 


In/Del 


2 A 


96 




GPNMB Glycoprotein NMB 


SNP 1 (BICF2P753624) [HpyCH4V] 


M1° 


26 


1.59 x 10" 4 


chrl 4:39,877,81 0-39,905,589 


SNP 2 (BICF2P1 34952) 


2 A 


1200 




MITF Microphthalmia transcription factor 


SNP 1 (BIGF2G630233682) [BspHI] 


2 A 


73 


3.80 X10" 5 


chr20:24,853,657-24,884,775 


SNP 2 (BICF2S23248988) 


1 A 


521 




TYR Tyrosinase 


Microsat 1 (GAAA) 1 7(GGAA)20(GAAA) 1 0 


2 § 


428 


1.00 x 10" 9 


chr21:1 3,797,070-1 3,891,31 7 


Microsat 2 (TTTQ10(TTTC)4(TTTC)13 


1 § 


70 




TYRP2 {DCT) Tyrosinase related protein 2 


SNP 1 (BICF2S23 137809) [sequenced] 


M2° 


33 


8.47 x 10" 4 


chr22:48,21 9,81 7-48,254,000 


SNP 2 (BICF2P452919) [Rsal] 


M1° 


8 




ASIP Agouti signaling protein 


SNP 1 (BIGF2P1 186810) [Mesl] 


1 A 


340 


3.57 x 10" 5 


chr24:26,327,360-26,366,307 


In/Del 


2 A 


105 




COMT Catechol-O-Methyltransferase 


SNP 1 (BICF2S22923369) [ApaLI] 


2 A 


82 


8.20 x 10" 5 


chr26:32,426,959-32,432,523 


Microsat 1 (TTTA)15 


1 § 


676 




GSK3B Glycogen Synthase Kinase 3-Beta 


Microsat 1 (CTA^)14 


M2* 


146 


1.00 x 10" 5 


chr33:26,516,949-26,699,712 


Microsat 2 (TTTA)13 


M2* 


143 





Key: 

a. Genes are listed in chromosomal order as obtained from the May 2005 build of the canine reference genome (UCSC Genome Browser: http://genome.ucsc.edu/), 
with the encoded proteins provided under the gene names. 

b. Microsatel lite-based markers are shown with the repeat type and number of perfect repeats present in the canine reference genome. Repeat blocks are 
separated by one to several nucleotides that do not match the perfect repeat. SNPs are listed with Broad Institute CanFam 2.0 SNP designation (http://www. 
broadinstitute.org/mammals/dog). One marker for ASIP is an in/del of an undefined but variable nature. The restriction enzyme used for PCR-restriction enzyme 
method of genotyping SNPs is shown in square brackets. Note some were genotyped by sequencing. 

c. The location of each marker is given with respect to coding region of each gene as seen in the 2005 canine reference genome on UCSC Genome Browser. 
Designations are upstream (1) or downstream (2) of the gene start site, or within the gene upstream of the exact midpoint (M1) or to the downstream of the 
midpoint (M2). 

d. Distances are given from each marker to furthest end of the gene from that marker. 

e. Each P value is the probability (combined probability of the two markers) of falsely excluding the true causative gene (see Additional file 2: Supplementary 
Methods for an example of how this was calculated). The method used for each marker is indicated in the column showing marker position: § used microsatellite 
mutation rate, A used rate of recombination between marker and coding region, o used rate of recombination from marker to nearest end of the coding region 
(probability of the recombination WITHIN the gene), * Clark's method for haplotype analysis [10]. 
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The exclusion analysis 

Exclusion analysis was performed under several assump- 
tions that are considered to be robust for canine Mendelian 
disorders; (1) all affected dogs share the same causal muta- 
tion, which is identical by descent due to a founder event 
(very likely, based upon previous results for common 
Mendelian mutations in purebred dogs; see Discussion), 

(2) no recombination has occurred between the marker 
and the mutation since the time the unique mutational 
event occurred (markers are very close to the gene and the 
probability of a rare event is calculated and presented), and 

(3) no new mutations have occurred in either the marker 
or the gene (also included in the probability calculation for 
a false exclusion). It is important to understand that, for the 
purpose of exclusions, only affected dogs are needed. This 
is analogous to exclusion in paternity testing. If an accused 
male does not share marker alleles with the baby, he cannot 
be the father; knowing or not knowing the marker allele 
frequencies in the population has no impact on the exclu- 
sion. It is only necessary to have allele frequency informa- 
tion when the alleles in the offspring are consistent with 
those of the accused male in order to provide a probability 
statement on how likely this match is compared to other 
possible males in the population. We sought to exclude 
most candidate genes as quickly and affordably as possible 
because only one gene can underlie a Mendelian disorder 
due to a founder event and the rest of the candidate genes 
are "innocent" with respect to the disease under exami- 
nation. An advantage to this approach is the cost savings 
obtained because it is not necessary to genotype unaffected 
controls. Controls only need to be genotyped if a gene can- 
not be excluded, in order to obtain positive evidence of an 
association between a locus and the disease phenotype. 
Note also that penetrance is not an issue because all Cairn 
terriers genotyped are affected and must share the neces- 
sary allele. 

With a likely dominant mode of inheritance, markers 
within or close to the gene harboring the causal mutation 
(or a necessary locus) would have at least one shared allele 
between affected dogs. For each of the 22 markers geno- 
typed there was no shared single allele between the OM- 
affected dogs thus providing evidence that none of the 
candidate genes underlie OM (Additional file 1: Table SI 
and Table S2). We then calculated the probability of falsely 
excluding the candidate genes: the p values ranged from 
1.59 x 10" 4 to 1.0 x 10" 9 (Table 1 and Additional file 2), 
making it very unlikely that any one the 11 candidate 
genes tested contains the OM mutation. 

Discussion 

Linkage-based and sequence-based exclusion analyses 
have been used in many studies to eliminate candidate 
genes hypothesized to underlie inherited phenotypes in 
dogs (for example: [11-13]). With late-onset conditions 



such as ocular melanosis compiling pedigrees (DNA 
and phenotypes from multiple offspring and at least one 
parent) can be difficult because it may not be possible to 
locate the siblings or the parents of affected animals, or 
they may already be dead. Sequencing of candidate genes 
can be time consuming and since the commonest ap- 
proach is to sequence only the exons and flanking por- 
tions of the introns, the method could miss causal 
mutations in non-coding regions. In this study we used 
association-based exclusion analysis, a technique that re- 
quires only the use of animals affected with the condi- 
tion and can be used for rapid exclusion of candidate 
loci. Obviously the hope is to identify a candidate locus 
that cannot be excluded and would thus warrant further 
investigation as the potential site of the disease causing 
mutation. Using association-based exclusion analysis we 
were able to exclude with high probability 11 candidate 
genes for OM in Cairn terriers. 

We argue here that the assumptions (no locus or al- 
lelic heterogeneity, or recombination and/or new muta- 
tions, as listed in the Results section) that underlie this 
approach are appropriate. Giger and colleagues noted 
that, with very few exceptions, the causative mutations 
that have been identified for canine genetic diseases are 
breed specific or affect very closely related breeds [14]. 
The relative rarity of a Cairn terrier-ocular melanosis- 
like condition in any other breed of dog makes it very 
likely that this condition also resulted from a single 
mutation event. 

We have estimated the probability of a recombination 
event causing a false exclusion of the true locus based 
upon the average recombination rate in dogs (1 cM/Mb). 
It has been suggested most chromosomal recombination 
occurs in "hotspots" in the human genome and 50% of all 
recombination occurs in less than 10% of the sequence 
[15]. In addition, recombination hotspots are relatively 
rare within or close to genes [16]. Although the region of 
linkage disequilibrium becomes smaller with increasing 
numbers of meiotic events (generations) after the founder 
mutation occurred, the timing of the first descriptions of 
ocular melanosis in Cairn terriers and the increasing inci- 
dence of glaucoma in the breed [17], suggests that the mu- 
tant allele became relatively common in the breed within 
recent time (probably less than 10 generations ago [2] and 
unpublished results). Thus, by using two markers, the 
average recombination rate becomes a conservative esti- 
mate, because of the small probability of the occurrence of 
two hotspots within or near a gene that would lead to a 
false exclusion. 

Although it has been commented in several publica- 
tions that some canine MS have observable mutations, 
to our knowledge no actual tabulation of mutation rate 
data has been published [18,19]. However, mutation rate 
data are available for humans, with the highest rates 
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equivalent to about 1 mutation per 100 meioses [20]. 
We have used this highest human mutation rate as our 
estimation of the rate of mutation at our MS. Although 
it is possible that canine MS mutation rates could be 
greater than this estimate, rates in dogs as little as two- 
fold higher would likely have been previously detected 
and published. With the approach used here, two rare 
events (recombinations and/or marker mutations) would 
have to occur to falsely exclude the causative gene. We 
also note for exclusion analysis, recombinations and mu- 
tations at loci which do not contain the disease causing 
mutation do not alter the conclusions inferred from the 
analysis; the only concern is for violations of the under- 
lying assumptions at the causal locus, leading to a false 
exclusion. With this approach and a single, necessary 
gene being responsible for the disease, most candidate 
genes can be screened in a relatively rapid and cost- 
effective manner. In this particular case, it is highly un- 
likely that the causal gene is among these 11 candidates. 

Conclusions 

This study provides strong evidence to exclude 11 po- 
tential candidate genes from association with the Cairn 
terrier OM locus. Further studies are required to identify 
the gene mutation that is causal for this condition. 

Methods 

Animals 

DNA was extracted by standard techniques from blood 
samples from privately owned pet Cairn terriers that had 
been diagnosed with OM by a veterinary ophthalmolo- 
gist. The clinical information provided was reviewed by 
SPJ or JTB. All procedures were in compliance with the 
ARVO statement for the Use of Animals in Ophthalmic 
and Vision Research and covered by Michigan State 
University IACUC AUFs 05-08-076-00 or 05-11-106-00. 

Development of polymorphic markers for candidate 
genes 

Candidate genes were selected because they either had a 
known involvement in ocular pigmentary disturbances 
in other species or a role in melanocytes and melano- 
some development. 

Two polymorphic markers were developed for each 
gene (Table 1). Microsatellites (MS) were identified from 
the USCS Genome Browser (http://genome.ucsc.edu/) 
May 2005 build of the canine genome using the repeat 
masker track. Single nucleotide polymorphism (SNP) 
markers were selected from the list of canine SNPs 
published by the Broad Institute of the Massachusetts 
Institute of Technology (http://www.broadinstitute.org/ 
mammals /dog). With the exception of COMT MSI and 
GPNMB SNP2 all markers were within 600 kb from the 
furthest end of the candidate for which they were 



selected. Primers to amplify each marker by polymerase 
chain reaction were designed using the Primer3 program 
(http://frodo.wi.mit.edu/primer3/). A universal labeling 
method was employed to fluorescently label the MS 
PCR products [21]. The MS markers were sized on an 
ABI PRISM 3130 genetic analyzer (Additional file 1: 
Table SI). Restriction enzyme digest assays were 
designed to allow genotyping of several of the SNPs (re- 
striction enzymes shown in Table 1) and the other SNPs 
genotyped by direct sequencing of the PCR products 
(Additional file l:Table S2). 

Calculation of probability of falsely excluding a genetic locus 

Markers close to the causal gene mutation are likely to 
be identical by descent in all affected alleles unless a re- 
combination event between marker and mutation oc- 
curred or the marker mutated. Either event could lead 
to falsely excluding the candidate locus. Therefore, the 
probability of a false exclusion of the causal gene due to 
such inconsistencies was calculated (see Additional file 2 
for an example of calculations). We first counted the 
number of necessary recombination/mutation events to 
obtain the genotype that would result in a false exclusion 
of the true culprit gene. Then the greater of the two 
probabilities of a recombination and a marker mutation 
was used for subsequent calculations. In all cases, if a 
marker was a SNP or an insertion/deletion the probabil- 
ity of a recombination was used (i.e., a recurrent SNP 
and insertion/deletion mutation is assumed to be always 
rarer than the probability of a recombination event for 
the data reported here). The probability of a recombi- 
nation event was assumed to be proportional to the dis- 
tance of the marker to the end of the most distant exon. 
The recombination rate per unit of DNA is equivalent to 
1 cM per Mb in the dog, based upon the estimates that 
the autosomal sex-averaged linkage map is 26.5 Morgans 
[22], and the physical map is 2631 Mb [23]. For all MS, 
the probability of a marker mutation was assumed to be 
0.01, based upon the mutation rate for a canine MS set 
at the highest mutation rate seen among human MS 
[20]. The probabilities for each recombination or marker 
mutation event were then multiplied for both markers 
if they flanked the gene on opposite sides, or for the 
number of steps for a single marker if it was contained 
within the gene (note that for markers that do not flank 
both sides of a gene, it is possible for a single recom- 
bination to produce non-shared alleles for both markers, 
although the much tighter linkage for markers contained 
within a gene can often make up for this disadvantage). 
The probability was then adjusted for the estimated 
number of generations (independent chances for a re- 
combination/mutation) since the foundation event; 1 - 
(1 - p) n , where n equals the estimated number of 
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generations since the foundation event (in this case, ten 
generations - based on first descriptions of OM). 

Additional files 



Additional file 1: Table SI. Genotyping results for microsatellites on 
the ABI 3130 genetic analyzer. Dogs with one recombination or mutation 
event are boxed. Light gray boxes indicate a second recombination or 
mutation event. Dark gray boxes indicate a third recombination or 
mutation event. Table S2. Bia I lei ic genotyping results. Bia I lei ic 
genotyping results with SNP information, SNP # correlates with Table 1, 
SNP position. 

Additional file 2: Supplementary Methods. Example of p value 
calculations using microsatellite markers. 
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